Electronic structure of random binary alloys : an augmented space formulation in 

reciprocal space 
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We present here a reciprocal space formulation of the Augmented space recursion (ASR) which 
uses the lattice translation symmetry in the full augmented space to produce configuration averaged 
quantities, such as spectral functions and complex band structures. Since the real space part is taken 
into account exactly and there is no truncation of this in the recursion, the results are more accurate 
than recursions in real space. We have also described the Brillouin zone integration procedure to 
obtain the configuration averaged density of states. We apply the technique to NisoPtso alloy in 
conjunction with the tight-binding linearized muffin-tin orbital basis. These developments in the 
theoretical basis were necessitated by our future application to obtain optical conductivity in random 
systems. 

PACS numbers: 71.23.-k 



I. INTRODUCTION 

The augmented space recursion carried out in a minimal 
basis set representation of the tight-binding linearized 
muffin-tin orbitals method (TB-LMTO-ASR) has been 
proposed earlier by us 0,0 as an interesting technique 
for incorporation of the effects of configuration fluctu- 
ations about the mean-field (the coherent potential ap- 
proximation or the CPA) for random substitutional^ dis- 
ordered alloys. This can be achieved without the usual 
problems of violation of the Herglotz analytic proper- 
tics 3] of the approximated configuration averaged Green 
functions for the Schrodinger equation for these random 
alloys. Earlier we had used this technique to look at 
short-ranged ordering in such systems 0, |5j , as well as 
local lattice distortions caused by size difference between 
the constituents of the alloy |(|. 

One of the dissatisfying features of the method, and 
this has to with the recursion part, is the truncation of 
the continued fraction expansion of the Green function. 
Truncation in the configuration space part of the problem 
can be handled easily. We truncate out only those config- 
urations which occur with low probability and contribute 
to the tail of the continued fraction, ft is the trunca- 
tion in real space on which we do not have a controllable 
handle. Any truncation in real space means that our re- 
cursion has been carried out on a finite cluster and edge 
effects become important. Quantities which converge fast 
are integrals of the density of states multiplied by well- 
behaved functions of energy. We can also estimate the 
errors committed by truncating at a particular step 0. 
However, the errors in the density of states itself cannot 
be controlled. This is because even a small perturbation 



(like truncation after a large number of recursive steps) 
has a profound effect on the spectrum of the Hamilto- 
nian (see Haydock 24]). The problem of truncation has 
always been laid at the door of the recursion method. 

Is it not possible to modify the TB-LMTO-ASR in such 
a way that the truncation is carried out only in config- 
uration space ? One way of reducing the gigantic rank 
of the Hamiltonian in a real-space-labelled basis, is to 
go over to reciprocal space. In the k labelled basis, for 
a basis involving only s, p and d states, the operators 
in reciprocal space have rank 9. However, to do this we 
require lattice translational symmetry. In a random bi- 
nary alloy, for instance, this is not immediately possible. 
However, the full augmented space, which is the direct 
product of the real space spanned by the site labelled 
basis {R^} and the configuration space spanned by the 
configurations of the system, possesses translational as 
well as point group symmetries Q. Configurations of 
a binary alloy can be labelled by a binary sequence of 
and 1 (or f and J. if Ising models appeal to you more) and 
uniquely described by the cardinality sequence {C}, i.e. 
the sequence of positions where we have a I or a J, state. 
We had shown earlier that in the subspace spanned by 
the reference states {0}, in which the configuration aver- 
age is described, we have full lattice translation symme- 
try provided the disorder is homogeneous The same 
statement would be true if there is short-ranged order or 
local lattice distortions, provided the short-ranged order 
or local lattice distortions are probabilistically identical 
anywhere in the system. A consequence of this is that 
probability densities are independent of the site label and 
the configuration averaged quantity : 



J2 J2 cxp {*( k ' Ri ~ k ' ' R ^} < G(Ri,Rj,z) > = G(k, z) <5(k - k') 

Ri Rj 
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Based on this, we had proposed a TB-LMTO-Recursion 
in the reciprocal augmented space jlCj . The recursion, 
as we shall show subsequently, is entirely in configura- 
tion space for each k label. The truncation is also in 
configuration space alone and leads to calculation of the 
configuration averaged spectral densities. These spectral 
densities are not a bunch of delta functions, as in the 
case of ordered systems, but the complex self-energies, 



J 



in general both energy and k dependent, shift the peaks 
as well as and broaden them : leading to fuzzy, complex 
band structures. 

Although our method allows us to carry our augmented 
space recursion in reciprocal space, for many physical 
problems we need to carry our integration over the Bril- 
luoin zone. For instance, to obtain the density of states 
or optical conductivity [Tl| : 



« n(E) » = / « A(k, E) » 

Jbz 8tt j 



— Tr [j e/ f (k, E, u) < A(k, E) > J e/ / (k, E, w) f < A(k, E + u) > 



r 



Another contribution of this paper is to modify the tetra- 
hedral method of Brilluoin zone integration, so that we 
may carry out a smimilar integration technique for inte- 
grands which are smoother than the highly singular spec- 
tral functions of the ordered systems. The proposed Bril- 
luoin zone integration is closely related to that of Jepsen 
[l2l | or Lehmann for ordered systems. 



II. AUGMENTED SPACE RECURSION IN 
K-SPACE 



The augmented space recursion based on the tight- 
binding linearized muffin-tin orbitals method (TB- 
LMTO- ASR) has been described thoroughly in a series 

of articles [11 [3 H H IH HI . We sha11 introduce the 
salient features of the ASR which will be required by us 
in our subsequent discussions. 

We shall start from a first principles tight-binding lin- 
earized muffin-tin orbitals (TB-LMTO) |H H3 in the 
most- localized representation ((3 representation). This is 
necessary, because the subsequent recursion requires a 
sparse representation of the Hamiltonian. In this j3 rep- 
resentation, the second order alloy Hamiltonian is given 
by, 



H (2) = E„ + h - hoh 



where, 



h =J2 (C« - E v r) Vr + Y, E A R 2 S ««' A «' 2 T ««' 

R R R' 

o=^o fl P s (1) 
R 

C R , E„# , A R and o R are diagonal matrices in angular 



momentum space : 

C R = Crl Sll' , E„t? = E vR l Sll' 
A R = A RL 5 LL > , = o RL 5 LL , 

and Srri is a matrix of rank L max . Vr = and 
Trri = \R) (R' I are projection and transfer operators in 
the Hilbert space Ti spanned by the tight-binding basis 
Here, R refers to the position of atoms in the solid 
and L is a composite label {£, to, to s } for the angular mo- 
mentum quantum numbers. C, A and o's are potential 
parameters of the TB-LMTO method; o _1 has dimen- 
sion of energy and E„'s are the energy windows about 
which the muffin-tin orbitals are linearized. 
For a disordered binary alloy we may write : 



Crl = Gt n R + Cf (1 - n R ) 



A]( 2 L = (Ai) 1/2 n R +(A*y"(l-n R ) 

R L 



J/ 2 



of n R + of (1 - n R ) 



(2) 



Here {n R } are the random site-occupation variables 
which take values 1 and depending upon whether the 
muffin-tin labelled by R is occupied by A or B-type of 
atom. The atom sitting at {R} can either be of type 
A (n R — 1) with probability x or B (n R = 0) with prob- 
ability y. The augmented space formalism (ASF) now 
introduces the space of configurations of the set of bi- 
nary random variables {n R } : $. 

In the absence of short-ranged order, each random vari- 
able n R has associated with it an operator M. R whose 
spectral density is its probability density : 



p{n R ) = xS(n R - 1) + yS(n R ) 
= -- lira 3?m(ti?, | {(n R 

7T 5^0 



iS)I-M R y I ti?.) 

(3) 
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where M# is an operator whose eigenvalues 1,0 cor- 
respond to the observed values of n R and whose cor- 
responding eigenvectors {|1_r), \0r)} span a configura- 
tion space (f>R of rank 2. We may change the basis to 
{| 1r), I Ir)} (see appendix A) and in these new basis 
the operator M.r is : 

n R -> M R = xV\ + yV R + ^xy (T fl u + ) (4) 

Two new vectors span the space <p R . The full configura- 
tion space $ = Yl% 4>R is then spanned by vectors of the 
form | HIT I • ■ •)• These configurations may be labelled 
by the sequence of sites {C} at which we have a J.. For ex- 
ample, for the state just quoted {C} = |{3, 5, . . .}). This 
sequence is called the cardinality sequence. If we define 
the configuration | |f . . . | . . .) as the reference config- 
uration, then the cardinality sequence of the reference 
configuration is the null sequence {0}. 
The augmented space theorem ^(J states that 

« A({n R }) > = < {0}|A|{0} > (5) 

where, 

A({M R }) = J... J A({X R }) l[dP(X R ) 

P(A_r) is the spectral density of the self-adjoint operator 
Mr. 



Applying this to JSJ to the Green function we get : 



<0G(k,z)>= (k® {0}|(zI-H< 2 )) 1 |k®{0}) (6) 

where G and H^ 2 ) are operators which are matrices in 
angular momentum space, and the augmented k- space 
basis |k, L eg) {0}) has the form 



(1/y/N) expHk -R)\R® {0}) 

R 

The augmented space Hamiltonian H^ 2 ) is constructed 
from the TB-LMTO Hamiltonian 0H^ by replacing 
each random variable n R by operators M. R . It is an oper- 
ator in the augmented space W = TL (E> < I > - The ASF maps 
a disordered Hamiltonian described in a Hilbert space 
Ti. onto an ordered Hamiltonian in an enlarged space 
where the space "J is constructed as the outer product 
of the space Ti. and configuration space $ of the random 
variables of the disordered Hamiltonian. The configu- 
ration space $ is of rank 2 N if there are N muffin-tin 
spheres in the system. Another way of looking at H'- 2 *' is 
to note that it is the collection of all possible Hamiltoni- 
ans for all possible configurations of the system. 

The resolvent of the Hamiltonian can be expressed in the 
following way : 



(zl - H< 2 )) = (zl-C- A 1 / 2 S A 1 ' 2 + hoh) ' 



= A" 1 / 2 



zl- C 



C-E b 



S (A^oA 1 / 2 



C-E„ 



A" 1 / 2 



Expressions in bold are matrices in angular momentum 
space and other than S and H^ 2 ) and G all others are 
diagonal matrices. 

In the above expression, since 



A" 1 / 2 = J2 A(A-!/ 2 ) V R ® I + B(A- 1 / 2 ) V R ® V 

R 



where for any diagonal (in angular momentum space) 
operator V : 



A(V) = A(V L ) 5 LL , A(V L ) =xVf + y V? 
B(V) = B(V L ) 6 LL , B(V L ) = {y-x) (V? - V L B ) 
F(V) - F(V L ) 5 LL , F(V L ) = ^ (Vf - VS) 

we get : 

A- 1 / 2 |k®{0}) = A(A- 1 / 2 )|k®{0})... 

... + F(A- 1 / 2 )|k®{i?}) = |1} 

The ket |1} is not normalized and we define the normal- 
ized ket as, |1) = [A(A _1 )] _1 / 2 |1}. Then we may rewrite 
© as 
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<G(k,z)> = (1| (zl- A + B + F-S 
• • • + (J + S) 5 (J + S 

where, 



A 
B 
F 



^{A^A^j/AfA- 1 )} V R ®T 

R 

^{B^zI-^A-^^A- 1 )} Vr< 

R 

^{F^zI-^A-^/ACA- 1 )} Vr( 



and J = 3 a + J_b + Jf and = 0,4 + 03+0^ where 



In case there is no off-diagonal disorder due to local lat- 
tice distortion because of size mismatch : 



S = 0A(0^R 1 r 1/2 0Srr> 0A(0A r })- 1 / 2 T rw ®1 

R W 

This equation is now exactly in the form in which recur- 
sion method may be applied. At this point we note that 
the above expression for the averaged Gll{^-, z) is exact. 
The recursion method addresses inversions of infinite ma- 
trices ??. Once a sparse representation of an operator in 
Hilbert space, fv- 2 \ is known in a countable basis, the re- 
cursion method obtains an alternative basis in which the 
operator becomes tridiagonal. This basis and the rep- 
resentations of the operator in it are found recursively 
^hrough a three-term recurrence relation : 



K +1 } - H( 2 V„} - a„(k)K} - #(k)K_i} (9) 



Ja = 
Jb = 

J F = 

OA = 

ob = 
o F = 



J2{0A(0C-0E„)0A- 1 )/0A(0A~ 1 )} V R ( 

R 

{0B((0C - 0E„)0A- 1 )/0A(0A- 1 )} V R < 

R 

J2 {0F({0C- 0E v )0A' 1 )/0A{0A- 1 )} V R ® [T^ + T fl iT } 

R 

^{0A(0 O 0A) 0A(0A" 1 )} V R ®1 

R 

^{0B(0 O 0A) 0A(0A~ 1 )} V R ®V R 

R 

^{0F(0 O 0A) 0A(0A~ 1 )} V R ® \jl 



with the initial choice |ui} = \RL) and 01 = 1. The 
recursion coefficients a n and (3 n are real and are obtained 
by imposing the ortho-normalizability condition of the 
new basis set as : 



a„(k) 



{n|H< 2 )|n} 



{n- l|H( 2 )|n} 



-[ 1 



T, 



{n\n} ' ; {n\n} 

and also {mlH*- 2 -* \n} = for m ^ n, n ± 1 

To obtain the spectral function we first write the configu- 
ration averaged L-projected Green functions as continued 
fractijSds : 



< G LL (k,z) > = 



- aii(k) - 



a2i(k) 



z - a NL (k) - r L (k, z) 



where T l (k, z) is the asymptotic part of the continued 
fraction. The approximation involved has to do with 
the termination of this continued fraction. The coeffi- 
cients are calculated exactly up to a finite number of 
steps {a n , P n } for n < N and the asymptotic part of the 
continued fraction is obtained from the initial set of co- 
efficients using the idea of Beer and Pettifor terminator 
[23 • Haydock and coworkers [2j| have carried out exten- 
sive studies of the errors involved and precise estimates 
are available in the literature. Haydock (2(| has shown 
that if we carry out recursion exactly up to N steps, the 



resulting continued fraction maintains the first 27V mo- 
ments of the exact result. 



It is important to note that the operators A, B,F are 
all projection operators in real space (i.e unit operators 
in k- space) and acts on an augmented space basis only 
to change the configuration part (i.e. the cardinality se- 
quence {C} ). 
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A\\{C}) = Ai||{C}> 

B||{C}> = A 2 \\{C})S(Re{C}) 

F\\{C}) = A 3 \\{C±R}) 

The coefficients A\ — A3 can been expressed from equa- 
tion (J7J. Similar expressions hold for the operators in 
equation JSJ. The remaining operator S is diagonal in 
k- space and acts on an augmented space only to change 
the configuration part : 

S||{C}) = £expHk. X )||{C-x}> 

x 

Here x s are the near neighbour vectors. The operation 
of the effective Hamiltonian is thus entirely in the con- 
figuration space and the calculation does not involve the 
space 7i at all. This is an enormous simplification over 
the standard augmented space recursion described ear- 
lier 0, 0, , where the entire reduced real space 
part as well as the configuration part was involved in the 
recursion process. Earlier we had to resort to symmetry 
reduction of this enormous space in order to make the 
recursion tractable. Here the rank of only the configura- 
tion space is much smaller and we may further reduce it 
by using the local symmetries of the configuration space, 
as described in our earlier letter [l^. However, this ad- 
vantage is offset by the fact that the effective Hamilto- 
nian is energy dependent. This means that to obtain the 
Green functions we have to carry out the recursion for 
each energy point. This process is simplified by carrying 
out recursion over a suitably chosen set of seed energies 
and interpolating the values of the coefficients across the 
band. 



where, 

< A kL {E) > = — lira {3™ < G LL (k,E-iS) >} 

7T (5^0+ 

To obtain the complex bands for the alloy we fix a value 
for k and solve for : 

z-E L {k)-X L {k,E) =0 

The real part of the roots will give the position of the 
bands, while the imaginary part of roots will be propor- 
tional to the lifetime. Since the alloy is random, the 
bands always have finite lifetimes and are fuzzy. 



IV. INTEGRATION IN K- SPACE 

To obtain the density of states we need to integrate over 
the Brillouin zone 

«n(£)» = £jf^_ (10) 

For ordered systems the spectral function is a bunch of 
delta functions : A^(E) = J2j A jK E ~ • B j( k ))> with 
j labeling a particular energy band. The integrand be- 
ing highly singular, the i nteg ral (|1(J|) has to be calculated 
carefully. Jepsen et. al. [12j and Lehmann 13] had pro- 
posed an accurate technique : the tetrahedron method, 
for obtaining such integrals accurately. In this section we 
shall discuss an extension of that method for application 
to disordered systems. 

In the presence of disorder the spectral function is 
smoother and we may rewrite it in terms of the real and 
imaginary parts of the disorder induced self-energy : 



III. SPECTRAL DENSITY AND BAND 
ENERGY 

The self-energy which arises because of scattering by the 
random potential fluctuations is of the form : 



/3 2 2 Jk) 



"2i(k) 



a NL (k)-T L (k,z) 



So the continued fraction can be written in the form 
l/(z-£ L (k)-£ L (k,£)), where E L (k) = a 1L (k). 
The average spectral function <C A^(E) 3> is related to 
the averaged Green function in reciprocal space as : 



« A k (E) > = < A kL (E) » 



E-E L (k)-ZX(k,E)) +££(k,£)2 



< A UL (E) > 



(11) 

Such a function is peaked around the zeroes of E — 
E L (k) - Ef(k,E) and the S£(k,E) provides the width 
of the peaks. The spectral function behaves roughly as 
Lorentzian in the vicinity of its peaks. We may reduce 
the above expression to one amenable to the tetrahedron 
integration form by the following trick : 



dE> Zgii^ 

(E-E>-ZK(k,E)) 2 + Z{(k,Er 

= I dE' W kL {E,E') 6 (E' - E L (k) 



s(e' 



Erik) 



where WkL is defined as a weight function. Now inte- 
grating above over the Brillouin zone, we may get con- 
figuration averaged density of states (DOS) : 



f <i 3 k f 

« n {E) > = ^ / dE ' ^ A kL(£) » 

At this stage, in order to simplify notation we shall drop 
the L index from all L dependent factors and understood 
that the eventual result is summed over all L. In order 
to perform the above integration over BZ, we have gen- 
eralized tetrahedron method developed by Jepsen et al 
[l2| and Lehmann et al 0] to include the weight func- 
tion W k (E, E'). We have followed the idea of MacDonald 
et al in this generalization the energies as well as 

the weight functions are linearly interpolated throughout 
the vertices of small tetrahedrons. We label the energies 
at the vertices of the ith tetrahedron E\,E\, E\ and E\, 
where the indices correspond to increasing magnitude of 
the energy, i.e., E\ > E\ > E\ > E\ and the corner val- 
ues of the weight function be W{ , W|, W| and W|. Then 
the averaged DOS may be written as : 

„ N 4 
« n(£) »= Fmz / dB' E CI E 7 * W * ( 12 ) 

i=l fc=l 

where P k = I k (E, E', E[ , E\, E\, E\), N is the number 
of tetrahedral micro-zones and Vmz is the micro-zone 
volume and also, 

for E\<E' < E\ 



c n = 


- 3 F2iF 3 iF i i/(E' — Ei) 


I\ - 


-- (F 12 + 


F 13 + F u ) /3 


II = 


-- Ffci/3, 


k = 2,3,4. 


E' < E\ 






& = 


3 (^23^31 + ^32^24) /E41 


I\ = 


Fu/3 + 


F13F31F23 / C l E±i 


P - 

1 2 — 


F23/3 + 


F^F^/C'Em 


P — 


F 32 /3 + 


^31-^23 / C l Fil 


P — 


F41/3 + 


F42 F '24F32 / C l E41 


E' < E\ 






C l = 


- 3 F14F24F34/ (E4, — E) 


P - 


= Ffc4/3, 


k = 1,2,3 




= (F41 + 


F 42 + F 43 )/3 



where E mn = E m - E n and F mn = (E' - E n )/E, 
Also < n(E) > is zero for E' < E\ or E' > E\. 




Number of shells 

FIG. 1: Showing how the number of lattice sites increase 
with increasing the number of shells in the real space and 
reciprocal space map. 



V. COMPUTATIONAL DETAILS AND 
RESULTS 

For ordered faces the calculations have been performed 
in the basis of linear muffin-tin orbitals in the atomic 
sphere approximation including combined corrections. 
The scalar-relativistic calculations in this case are car- 
ried out for equal atomic spheres. The k- space integra- 
tion was carried out with 16 x 16 x 16 mesh resulting 145 
k-points for cubic primitive structure in the irreducible 
part of the Brillouin zone. 

In Fig. ^ shows how the size of the augmented space 
map (in both k- space and real-space representation) in- 
creases as we increase the number of nearest neighbour 
shells from a starting site. We note that the reciprocal 
space map at a particular recursion step is much smaller 
than the real augmented space map. This is because 
in the reciprocal augmented space we generate only the 
different configurations. The full real space lattice map 
has been collapsed using lattice translational symmetry 
in full augmented space. 

We have first carried out calculations on a simple model 
disordered binary alloy system described by a s-state 
tight-binding Hamiltonian with nearest neighbour hop- 
ping integrals only. In Fig. [21 we compare the results 
obtained using reciprocal and real space formulation of 
ASR. The k- space integration has been performed in 
two ways. The brute force method, where we replace the 
integral by a sum with appropriate weights at different 
k-points, generates some unusual oscillations particularly 
in the lower part of the band. However, the tetrahedron 
method gives smoother results which are in good agree- 
ment with the real space calculations as well. 
We now go over to calculations for the disordered 
Ni5oPt5o alloy. We have used the minimal basis set of 
the TB-LMTO with nine orbitals per atom (s, p and d) 
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FIG. 2: Comparison of the average (50-50) density of states for a model fee alloy calculated using k- space formulation of 
ASR (solid curve) with the same calculated using real-space formulation of ASR ((dotted curve), k- space integration has been 
performed in two ways : (a) using Tetrahedron Method (right solid curve) (b) multiplying spectral function Ak(E) by k-point 
weight and then summing up over k (left solid curve). In both figures we note that the oscillations shown by the brute force 
technique is smoothened by the TM. 



to set up our Hamiltonian. 

In Fig. |31 we present the results for the spectral functions 
for Ni5oPt5n alloy along the Y — X direction. We have 
chosen 11 k-points having equidistant between r to X 
points and have shown spectral function in those points. 
These spectral functions shows good agreement with the 
same results obtained from KKR-CPA calculations [? ]. 
If we look carefully, we see that the widths of the spectral 
function varies considerably as a function of k and E. 
There are some simple trends concerning this behaviour. 
The sharp peaks on the lower band edge of near the T 
point appear as the s-like band. As we go towards Y to 
X point the s-band hybridizes with the p-band and the 
peak becomes wider. The structures on the upper band 
edges are mostly due to the overlap of the d-states of Ni 



i 50- 




-0.5 

Energy (Ryd.) 



FIG. 3: The spectral function of NisoPtso alloy plotted as a 
function of energy at several k-points along the F — X direc- 
tion. 



and Pt. Disorder affects to these d-dominated states are 
strong and there is significant broadening. 

In Fig. 0] we presented the complex fuzzy bands of the 
disordered alloy. The disorder smearing is maximum in 
the overlapping d-bands of the constituents, and is neg- 
ligible in the s-like part. This is also apparent in the 
spectral functions shown earlier. The sharp s-like peaks 
flank wide d-like structures in Fig. [3] 
Finally using our modified tetrahedron method we have 
calculated the density of states (DOS) of ordered and ho- 
mogeneous disordered NiPt alloys from its spectral func- 
tion. Side by side we have also carried out the same 
calculation in real augmented space. In Fig. [3] we show 
the I— projected density of states for the Ni 50 Pt5o alloy. 
We compare the k-space results with those found from 
real space recursion. The main improvement occurs in 
the eg and t2g d-bands. In particular, the sharp feature 
straddling the Fermi energy is better reproduced in the 
k-space recursion than that in real space. The reason 
for this is the early truncation of recursion in real space 
and the consequent finite size effects to which the more 
localized d-states are more susceptible. 

In Fig. (top row) we show a comparison between aver- 
age DOS calculated by real and reciprocal-space recur- 
sions. As discussed before, it is the sharp feature strad- 
dling the Fermi energy with major contribution coming 
from the NMi d-states which are not well reproduced in 
the real space technique. In this point our k- space calcu- 
lations agree with the KKR-CPA results of Staunton et. 
al. 13 • In the left lower panel of Fig. Elwe show the DOS 
for pure Ni and Pt, but in a lattice with the lattice pa- 
rameter same as the alloy. We may compare this with the 
DOS for the disordered alloy. The right most three peaks 
at -0.25, -0.16 and -0.11 Ryd. of the disordered DOS are 
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FIG. 4: (Right) Ni and Pt energy bands on a lattice appropriate to the NisoPtso alloy, in the V — X direction. Average lattice 
parameter ao = 7.127ou was fixed after minimizing the energy. (Left) The fuzzy band of the disordered NisoPtso system plotted 
along the same direction. 



mostly contributed by Ni whereas the left (lower energy) 
structures come (large shoulder at -0.57 Ryd.) mostly 
from Pt. The sharp peaks in the elemental results are 
obviously because of the Van Hove singularities of the 
DOS. The effect of disorder mainly smears out the sharp 
peaks present in the DOS . The disorder smearing is more 
pronounced for the d-like parts of the band. We remark 
that there is very little shift in the DOS-related features 
between the ordered and disordered states. 
Finally in the right lower panel we show the photoem- 
mission spectrum of NisoPtso reported by Nahm et.al. 
|28| . The general features with a double peak straddling 
the Fermi energy and a lower energy shoulder are clearly 
seen. The photoemmission spectra are convolutions of 
the density of states with a weakly energy/wavenumbcr 
dependent transition matrix. This may lead to shifting 
and smearing of the prominent peak structures. Keeping 
this in mind, our k-space recursion results are in good 
agreement with experiment. 

VI. REMARKS AND CONCLUSION 

We have presented here an augmented space recursion 
formulation in reciprocal space. We also present a gener- 
alization of the tetrahedron method proposed by Jepsen 
et. al. 12] for inverting the spectral functions to obtain 
the density of states. This technique will be useful for 
carrying of Brilluoin zone integrals for disordered alloys. 
We have studied both a model alloy and NiPt. The latter 



was so chosen as its has a sharp structure straddling the 
Fermi energy and will be a sensitive test for the accuracy 
of our technique. 
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APPENDIX A: THE AUGMENTED SPACE 
THEOREM 

Let f(n R ) be a function of a random variable n R , whose 
binary probability density is given by : 

p(n R ) = x 5(n R -l)+y S(n R ) 

We may then write : 

p(n R ) =-- lim 5Sm(\ R | ((n R + iS)I - Mr)' 1 \ U) 

Here, the operator Ts/L R acts on a space spanned by 
the eigenvectors and \0 R ) of M#, corresponding 

to eigenvalues 1 and 0; | T-r) = v^llfi) + y/y\Qn) is 
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FIG. 5: Comparison of the partial density of states of NisoPtso alloy calculated using augmented space recursion in (a) 
real-space formulation (left panel), (b) k- space formulation (right panel). 



called the reference state. Its orthogonal counterpart is 
I Ir) = \/v\^r) — Vx\0r). The representation of M. R in 
this new basis : 



M 



R = 



X 

xy 



/xy 

y 



Now, 



< /(»«) > 



f 



f(n R )p(n R )dn R 



f(n R )(U |(n fl I-Mfl) x | U) dn R 



= --^ m E E /(««)<T«|A>... 

A=0,1 A'=0,1 

...{AKnal-MflJ-^A'XA'l U) dn R 
= E tt* WfWiM Tfl) 

A=0,1 

= (Ttv. |f| T«> 



(Al) 



Here f is an operator built out of f(n R ) by simply re- 
placing the variable n R by the associated operator M. R . 
The above expression shows that the average is obtained 
by taking the matrix element of this operator between 
the reference state | The full Augmented Space 

Theorem is a generalization of this for functions of many 
independent random variables {n R }. 



APPENDIX B: TERMINATORS 

The recursive calculation described earlier gives rise to 
a set of continued fraction coefficients {a„,/3„}. In any 
practical calculation we can go only upto a finite num- 
ber of steps, consistent with our computational process. 
In case the coefficients converge, i.e. if \a n — a\ < e, 
\0n — j3\ < e ior n > N 7 we may replace {a n ,(3 n } by 
{a, 0} for all n > N. In that case the asymptotic part 
of the continued fraction may be analytically summed to 
obtain : 
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FIG. 6: Comparison of the density of states NisoPtso alloy calculated using augmented space recursion in (a) k- space 
formulation, (b) real space formulation, (c) Density of states of Ni (solid line) and Pt (dotted line) on a lattice appropriate to 
the NisoPtso alloy, (d) Valence-band photoemission spectra of NisoPtso with photon energy hv = 60. 



T(E) = U E -a- ^(E-af-^A 

which gives a continuous spectrum a — 2(3 > E > a + 2(3. 
Since the terminator coefficients are related to the band 
edges and widths, a sensible criterion for the choice of 
these asymptotic coefficients is necessary, so as not to give 
arise to s pur ious structures in our calculations. Beer and 
Pettifor [23 suggest a sensible criterion : given a finite 



number of coefficients, we must choose {«,/?} in such a 
way so as to give, for this set of coefficients, the minimum 
band width consistent with no loss of spectral weight 
from the band. Let us call these values {a c , (3 C }. This cri- 
terion is easily translated into mathematical terms. The 
delta function that would carry weight out the band must 
then be situated exactly at the band edges. We thus de- 
mand that the continued fraction diverge simultaneously 
at both the top and the bottom of the band. 

At the band edges : T(a ± 2(3) — ±/3 and so, 



<G(a±2/3)>= 



A?/2 



#/4 



±(3 — -(«2 — ") 



±/3 — {oln — ot) 



For a given a, the (N 
agonal matrix 



1) eigenvalues of the finite tridi- 



/ |(ai - a) ±/? 2 
\(3 2 §(a 2 -a) I As 



\ 
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are values at which the Green function diverges. The 
maximum and minimum of this set of eigenvalues are 
those values of (3 for which spectral weight has just split 
off from the band. Thus our choice of a is that value 
for which the maximum eigenvalue is the largest and the 
minimum the smallest. Since the terminator only in- 
volves [3 2 we must have 



/3 C = sup I3 max {ac) = inf |/3 m i„(a c )| 

{a} W 

With this choice the terminator T(E) has all the Herglotz 
properties required. 
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